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ABSTRACT 

The relationship between longitudinal development of 
muons and conventional equi-intensity cuts is carefully 
investigated. The development of muons in EAS has been 
calculated using simulation with a scaling violation 
model at the highest energies and mixed primary 
composition. Profiles of equi-intensity cuts expected 
at observation altitudes of 550, 690 and 930 gem*’ 2 
can fit the observed data very well. 

1. Introduction In recent years estimates of the longitudinal 
development of muons in EAS have been presented by several groups, based 
on equi-intensity cuts of muon size spectra following the method 
introduced 20 years ago by the BASJE group [1]. As is well known the 
relationship between such equi-intensity cuts and the profile of 
longitudinal development of the muons is not as straightforward as in 
the case of the electromagnetic component, and is particular to the 
altitude of observation. This arises because the muon longitudinal 
development itself depends on the zenith angle and the nature of the 
primary particle, moreover observations are made in general above a muon 
threshold energy which is zenith angle dependent. Some reported 
calculations when compared with measurements appear to show significant 
discrepancies [2], [3]. We report here calculations of equi-intensity 
curves, based on simulated showers, for different fractions of Fe nuclei 
in an otherwise pure proton beam using a model for hadronic interactions 
consistent with the degree of scaling violation discussed in [4]. 

2. Model Used in the Simulations A model for p-p collisions was 

adopted which incorporates scaling violation for interactions of the 
leading particle at E > 2 10 13 eV to the extent proposed in [41. with 
radial scaling assumed for all other interactions. Above lCr^eV the 
cross-section for hadron air-nucleus interactions was assumed to rise as 
0 mel. = (l+ a £n 2 E) , with values of the interaction lengths below this 

energy given by A air = 90 gem' 2 and A^^ = A k _ ai =120 gem -2 . For 
heavy nuclei incident we assumed u^-air = TTr o^ 1/3+ air • Th e 
smaller elasticity [5] and larger multiplicity [6] reported for p-air 
nucleus collision was allowed for by taking P p _ a i r = 0.31 and 
<m p-air > = l*^^ <in pp > - Although these additional particles prising from 
a nuclear target are in the central region and contribute little to the 
shower size they make a noticable contribution to the muon component at 
low energies. A T distribution was taken for the transverse momentum, 
whose mean, <p t >, was assumed to vary with the density of particles in 
rapidity space, consistent with recent accelerator data [7]. 

Three dimensional Monte Garlo simulations were made for proton 
and Fe nucleus primaries at different energies, at the vertical and at 
zenith angles of 30°, 40°, 50° and 60°. In fig.l we show the average 



HE 4.3-12 


i3e 



Expected muon 
longitudinal ; . 
development at 
fixed primary 
energies initiated 
by protons (solid 
curves) and iron 
dashed curves) . 

The dot-dashed 
curve is for radial 
scaling model. 


vertical longitudinal development of muons, N y (E o ,0=O,x,e > 0.6GeV) at 
fixed energies, the case for protons where radial scaling is assumed at 
all energies is also shown for comparison; e here is the muon threshold 
energy. These curves are generally similar to the calculations reported 
in [3]. 

3. Calculation of Equi-intensity Cuts To replicate the method used 
to derive equi-intensity cuts in the experiments would require the 
derivation of muon size spectra at different depths and zenith angles 
arising from particles selected from a primary spectrum of given slope 
and composition. This would involve an inordinate amount of computation, 
instead an approximate method is used. Because the muon decay proba- 
bility depends on zenith angle, vertical longitudinal profiles cannot 
be used for constructing equi-intensity curves for comparison with 
experiment, unlike the case for the electron component. Longitudinal 
developments at different zenith angles at each depth^must be used, 
in addition, as already mentioned, the relevant muon threshold energy 
must, in general, be taken as e o Sec0 where e 0 is the threshold for 
vertical muons. At a vertical depth Xq, when there is only one primary 
species present the experimental equi-intensity cuts correspond to 
N y(E o ,0,x = x o Sec0,e>e o Sec0) — assuming a unique relationship between 
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tne energy of the primary E 0 and muon size at 0, J(Ny,x) = J(E 0 (N u ,x>) 
(fluctuations may be allowed for by considering instead root mean square 
sizes [8]). If we consider a primary beam containing two species, protons 
and Fe nuclei say, we have 

J ( N y»x) = J p (E^ p) (N^,x)) + J Fe (E^ Fe) (N y ,x)). 

If the primary spectrum J(E 0 ) = AE Q ^ is assumed to contain protons, with 
a constant fraction k, and Fe nuclei only an equi-intensity curve is 
defined by 

AE^-Ytk+d-kJCE^/E^hh = Jy, a constant...,. (1) 


Since the range of primary energies of one species contributing to any 
one. equi-intensity curve will not be very great we approximate 
Ny (E 0 ,X) = E^gp(x), Ny e -'(Eg,x) = Eogp e (x), where w may depend on 
depth or zenith Sngle. Defining A(x) = (g p (x)/g Fe (x)) 1/& and using (1) . 
^we may write Ny(x,J„) = rf 0 g p (x), where E£*= eJp>(N u ,x) [k+(l-k) A" Y rx ) ] l ' y 
E o (Ny>x) - (A/Jy) 'i. Thus when k is not 0 or 1 equi-intensity curves > 
should be constructed from the simulated longitudinal developments, 
at a fixed E Q but at values of E p which depend on the depth. 


not 


4 - Results Using the profiles simulated at different zenith angles 

values of the function A(x) could be obtained. Taking A = 2.5 x 1(T, 



Expected muon long- 
itudinal development for 
an integral intensity 
J = 2.5 x 10" 8 m" 2 s“ 1 sr” 1 
at Chacaltaya, Tien 
Shan and Akeno. The 
dot and dot -dashed 
curves correspond 
to 50% and 25% of 
iron in primary cosmic 
ray. 
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y = 2 [8], expected equi-intensity 
curves for Jy = 2. 5 10” 5 6 7 8 m'* 2 s” 1 sr' 1 at 
Chacaltaya, Tien Shan and Akeno have 
been calculated and are shown in fig. 2, 
where a significant difference from the 
profiles at constant energy can be 
observed. Curves for 75%p, 25%Fe, 
normalised in each case to give a best 
fit to the data, are compared to 
experimental data from these experiments 
in fig. 3. The relevance of these curves 
to the primary composition is noted in 
another paper (OG 5.2-12) at this 
conference. 
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